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INTRODUCTION 


The  Sensor  Accuracy  Check  Site  (SACS)  is  located  at  the  Long  Beach  Naval  Support 
Activity,  Long  Beach,  California.  The  facility  is  operated  by  the  Naval  Undersea  Warfare 
Engineering  Station  and  technical  direction  is  provided  by  Naval  Ocean  Systems  Center, 
SACS/FORACS  Group  (Code  47).  The  principal  purpose  of  SACS  is  to  measure  sonar 
transmitter  and  receiver  performance  for  Naval  combatant  ships.  The  ship  is  moored  with 
the  sonar  transducer  positioned  in  the  center  of  a  1 20-ft-radius  pier.  A  6-ft-long  test 
transducer  is  mounted  on  a  movable  carriage  which  can  traverse  the  perimeter  of  the 
circular  pier  for  collection  of  data. 

The  SACS  environment  includes  several  physical  phenomena  which  affect  performance 
measurements.  These  phenomena  must  in  part  be  studied  with  the  use  of  an  acoustic 
mathematical  model.  This  document  presents  initial  development  of  the  model  in  four 
sections.  Section  I  derives  equations  for  the  summation  of  two  independent  acoustic  waves 
at  a  point  receiver.  Section  II  derives  equations  for  the  summation  of  multiple  acoustic 
waves  at  a  point  receiver.  Section  III  defines  a  three-dimensional  source  cardioid.  Finally, 
section  IV  develops  general  equations  for  the  summation  at  a  point  receiver  of  direct  and 
reflected  path  acoustic  waves  emanating  from  multiple  cardioid  sources. 

This  is  an  intermediate  document.  A  more  extensive  technical  document  will  be 
published  at  a  future  date  as  the  model  is  further  developed. 
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SECTION  1.  SUMMATION  OF  DIRECT  AND  REFLECTED 
PRESSURE  WAVES  FROM  A  SINGLE  SOURCE 


TRIGONOMETRIC  DERIVATION 


This  section  provides  the  solution  for  summing  two  waves  at  a  nondirectional  point 
receiver.  The  first  wave  travels  a  direct  path  from  the  source.  The  second  travels  a  reflected 
path.  Both  waves  are  projected  simultaneously  from  a  directional  point  source  as  shown  in 
figure  1. 

If  the  reflected  path  length  from  the  source  to  the  receiver  is  rf,  the  direct  path 
length  from  the  source  to  the  receiver  is  rd,  the  sound  pressure  at  the  defined  reference 
distance  (eg,  1  yard)  in  the  direction  of  the  reflection  point  is  Ar,  and  the  sound  pressure 
in  the  direction  of  the  receiver  is  Ad,  then  the  forms  of  the  reflected  wave  and  the  direct 
wave  at  the  point  of  the  receiver  are  (see  reference  1 ) 


Ad  /  2jrrd  v 

"d — —  cos  ^cot - - — J  direct  pressure  wave 

\  /  2rrrj.  \ 

Pr  =  — ~ —  cos  ^cot - - —  J  reflected  pressure  wave 


where  k  =  reflection  coefficient  (minus  sign  indicates  a  180°  phase  shift  at  the  surface), 

X  =  wavelength, 
t  =  travel  time, 
w  =  angular  frequency. 

The  directionality  of  the  source  is  implied  by  allowing  the  coefficients  Ad  and  Ar  to 
be  unequal.  The  two  waves  are  traveling  toward  the  receiver.  This  is  observed  by  requiring 
the  quantity  (rot  -  2rrrd/X)  to  be  zero  which  corresponds  to  the  peak.  Then,  as  time  t 
becomes  larger,  rd  must  also  become  larger  to  keep  the  quantity  equal  to  zero.  Therefore, 
as  time  increases,  the  peak  moves  toward  larger  values  of  rd;  ie,  toward  the  receiver. 

s'i\ 

! 

I  ^  Q  REFLECTION  POINT 


DIRECTIONAL 

SOURCE 


SURFACE 


OMNI 

RECEIVER 


Figure  1 .  Interference  at  R  between  sound  originating  at  source  S  and  the  surface  image  S'. 


1.  Technology  of  Underwater  Sound,  p  13,  Report  NADC-WR-6509,  May  1965 


Interference  will  occur  at  the  receiver  since,  tor  example,  the  reflected  pressure  can 
vary  front  -Ar/rr  to  +Ar/rr  The  summation  of  the  two  pressure  waves  at  the  receiver  can 
be  defined  as 


The  intensity  is  defined  as 


I  = 


(pressure)- 

pc 


(5) 


where  p  =  medium  density  and  c  =  propagation  speed.  Since  only  the  relative  intensities  are 
of  interest,  the  factors  p  and  c  are  not  necessary  and  only  the  value  tor  the  pressure  square 
is  needed.  Therefore,  squaring  equation  (4).  we  have 


<Pd  +  pr>- 


Ad 
rd  rr 


2  W  ,  2ltr*\ 

rr  cos 


-  rdrrcos  (“'  -  TT  )  cos  (“‘  •  -T~) 

,k^r:cos2^t.i^) 


Using  the  trigonometric  identity 

cos  a  cos  P  =  ^  cos  (a  +  P)  +  cos  (a  -  (3) 


(6) 


(7) 


or 


cos 


/  2irrd\  /  2rrrr  \  ,  /  2ff(rd  +  rrM 

y«t-—  Jcos^wt-—  f  --^cos^.wt-  x  ; 

|  /  2ir(rr  -  rd)  \ 

+  -cos\ - x - )' 


(8) 
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and  substituting  in  equation  (6),  we  have 


Pr>2  *  -^2  ’’  C“2  (“’  ‘  T1)  -  V'  C“  (J“‘  ‘ 


2ir(rT  +  rd) 


A.  /  2ir(rr  -  rd)  \  -  A2  -  /  2trrr  \| 

-krdrdrrcos  (— x - )+k  ^  rd“»2  («*-—))•  (9) 

It  is  now  of  interest  to  determine  the  time  averaged  intensity  at  the  receiver.  Consider  the 
cosine  forms  in  equation  (9)  such  as: 

cos2  (cot  -  a) 


cos  (2cot  -  P). 

The  time  averages  of  these  are 

T 

—  J  cos2(cot  -  a)  dt 


t  =  time  of  one  full  period  . 


x  =  cot  -  a 
dx  =  codt 


dt  =  — dx 
co 


1  f  2  a  1  1  A  1  ■  , 

—  /  cos^xdx= —  —  x +—  sin  2x 
CJT  J  COT  2  4 

0 


=  —  •  -1  (cot  -  O)  +  \  Sin  2(cot  -  a) 
cot  2  4 


11  a  x  1  •  -v  ** 

=  —  -  cot  — ;  +  7  sin  2(cot  -  a) 

cot  2  2  4 

-  0  +  ^  --  sin  2(0  -  a) 

2  4 


Therefore, 


—  f  cos2  (cot -a)  dt  =  —  ^  cot  =  4  . 
r  J  cot  2  2 
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Next 


T 

—  I  cos  (2cot  -0)  dt  . 
T  0 

x  =  2cot  -0 
dx  =  2codt 

1 

dt  =  —  dx 
2  co 


r 


cos  x  dx  =  - — 
2  cor 


sin  x 


(12) 


1 

2  cor 


:cor 


sin  (2cot  - 0) 


sin  <  2cor  -  0)  -  sin  (0  -  0) 


(13) 


Therefore. 


—  J  cos  (2cot  -  0)  dt  =  zero. 

T  0 

Substituting  the  time  averages  of  equations  ( 1 1)  and  (14)  into  equation  (9),  we  have 


(14) 


<pd  +  Pr^ 


r3 


rr  ^r  (  27r^rr  ~  rd^  \  A"rd 

T-k^r«r'con  x  r-^r 


(IS) 


or. 


2  Ag  k2A2r  kAdAf  /  27r(rr-rd) 


(Pd  +  Pr)“  =  —  t - - 

2rd  2r2  rdrr 


■cos 


(  27Krr-rd)j 


(16) 


This  is  the  relative  average  intensity  at  the  receiver.  It  is  dependent  on  the  path  length  differ¬ 
ence  (rr  -  rd)  and  has  a  maximum  value  of  '/2(Ad/rd  +  kAr/rr)2,  which  occurs  when  rr  is  180° 
longer  or  shorter  in  phase  because  of  the  1 80°  phase  change  on  reflection  at  the  surface. 
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COMPLEX  NOTATION  DERIVATION 

The  summation  of  two  pressure  waves  may  also  be  derived  with  complex  notation. 
The  basic  wave  equations  are 


Ad  *(<^t 
*  rdC 


2trrrt  * 


(real  part 


(real  part 


direct  path) 


reflected  path)  (18) 


The  sum  of  the  pressure  waves  is  then 


+P  =2iei(“'  -^r) 

'  rd  rr 

./  27rrdv  /  2rrrrv 

+  P  =  ^e\"  X  re‘V  X  )  eitot 


rd  T  rT  ~\  — c  '  - - '  e- — .  uu) 

[  rd  rr  J 

The  time  average  intensity  of  the  sum  would  be  given  by 

(Pd  +  Pr)2  =-j(P<j  +  Pr)  •  (Pd  +  Pr)  *  (see  reference  2)  (21) 

The  bar  indicates  time  averaging .  The  asterisk  denotes  the  complex  conjugate. 

A  complex  number  of  the  form  Ae*^  can  be  thought  of  as  a  vector  in  complex  space  such  as 
shown  in  figure  2. 


i  IMAGINARY 
I  AXIS 
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REAL 

AXIS 


So  from  Filler’s  formula.  Ae‘y  =  A  cos  0  +  iA  sin  0  defines  the  two  components  of  the  vector. 
If  two  of  these  vectors  such  as  Pj  and  Pr  are  to  be  added,  then  from  equation  (20). 


or. 


(23) 


The  time  averaged  intensity  would  be  one-half  the  product  ot  equation  (23)  and  its  complex 
conjugate  as  defined  in  equation  (21 ).  (Note  e'401  •  e-'40*  =  e®,  which  becomes  1 .) 

Substitute 


cos  6  =  cos  (-0) 

(24) 

-sin  6  =  sin  (-0) 

(25) 

(26) 


(27) 


(28) 
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Then 


(pd  +  pr)24 

1 


(a  +  i/3)  •  (a  -  i/3) 
(a2  +  /32)  • 


(29) 


Substituting  (27)  and  (28)  into  (29)  and  expanding  the  squared  terms,  we  have 

/n".  n  ,2  1/Aa  ?(2md\  .kAdAr  /2,rr d\  t2*rr\ 

(Pd  +  Pr>~  =  2  V7C°S  \T')-2'^-COS\-X-/  COSlT-) 


rd 

k2A2 


“Ar  7/27rrr\\  lAX  /27rrr\ 

7coi~)  MTs,n  '~) 

*r  r 

kAdAr  '27rrd\  /2irTT\  2(2nTd\\ 

2-^7s,n(— ) s,n  l— )  *7  sm‘  (— ) ) 


Using  sin20  +  cos20  =  1,  we  have 


(30) 


- T  4  k2A2  kAdA  /  /2irrd\  /  2irr  \  /2irrd\  /  27rrr 

(P« +  r)  '4  +l^"^Tvos  W cosl— ) + sin  rr ) s,n  It1))  ■ 


(31) 


Now,  by  using  the  trigonometric  relation 

cos  x  cos  y  +  sin  x  sin  y  =  cos(x  -  y), 
we  obtain  the  following. 


(pd  +  Pr)2  = 


4 

2ri 


k2A2  kAdAr 

+  — - - cos 

2r2  rdrr 


^27r(rr-rd) 


)■ 


(32) 


This  is  the  same  as  the  value  for  the  relative  average  intensity  in  equation  ( 1 6)  found  by 
using  trigonometric  expressions.  The  complex  notation  solution  redefined  Pd  and  Pr  as  vectors 
illustrated  in  figure  3.  Referring  to  equations  (17)  and  (18),  the  vectors  have  magnitudes  of 
Ad/rd  and  -kAr/rr  The  angles  <f>  and  i p  are  the  angles  (coT  -  27rrd/\)  and  (toT  -  2nrr/\). 

The  two  vectors  are  then  added  as  shown  in  figure  4.  The  terms  in  brackets  in  equation  (23) 
are  the  sums  of  C  and  D  on  the  real  axis  and  the  sums  of  E  and  F  on  the  imaginary  axis. 

(Pd  +  Pr)2  in  equation  (29)  is  side  G2  plus  side  H2  which  is  (Pd  +  Pf)2  in  figure  4;  the  time 
averaging  is  accomplished  with  the  1/2. 
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>  IMAGINARY 


Figure  3.  Geometric  representation  of  complex  numbers. 


Figure  4.  Graphical  addition  of  complex  numbers. 


Therefore,  the  time  averaged  intensity  summation  of  a  direct  and  a  reflected  path  from 
one  source  is  expressed  as 

'  ^  kAHAr  /^trr-rH>\ 

(33) 

rHrr  \  a  / 
r 


- 7  Ag  k-A-  kAdAr  /2ir(rr-rd)^ 

(p.i  +  prr  =  — ;  +  — r  - cos  I  — - 1 

d  r  2rJ  2r2  rdrr  V  X  ' 


where  Ad  and  Ar  will  be  determined  by  the  cardioid  of  the  source  pressure  and  rd,  rr  will  be 
determined  by  the  geometry  of  the  problem. 
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RELATIVE  AVERAGE  INTENSITY,  dB 


An  example  of  the  relative  intensity  for  a  source  and  a  receiver  moved  in  depth  from 
the  surface  to  50  ft  is  shown  in  figure  5.  The  peaks  correspond  to  a  condition  at  which  the 
quantity  (rr  -  r^)  is  equal  to  a  multiple  of  one  wavelength  (1  ft).  The  peaks  are  broader  at 
lower  depths  because  the  rate  of  change  of  (rr  -  r^)  decreases  as  the  depth  becomes  larger. 

A  plot  of  (Pj  +  Pr)2  as  a  function  of  range  for  a  fixed  depth  is  shown  in  figure  6. 
Here  the  inverse  square  spreading  loss  is  more  easily  observed  and  again  the  rate  of  change  of 
(rr  -  r^)  decreases  as  the  range  increases. 

For  both  these  plots  a  cardioid  of  the  type  shown  below  was  used  for  the  source 
pressure.  The  MRA  (Maximum  Response  Axis)  of  this  cardioid  was  placed  in  a  plane 
parallel  to  the  surface  and  at  the  depth  of  the  receiver. 


cardioid  in  polar  coordinates  (p,  0)  (34) 

p  =  Vi  ( 1  +  cos  0 )  . 


x  -  i  o  tt 

k  -  1.0 


Figure  5.  Relative  average  intensity  vs  depth  for  a  single  source  and  a  point  receiver 
at  a  fixed  distance. 
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RELATIVE  AVERAGE  INTENSITY,  dB 


♦« 


-6 


SOURCEDEPTH  -  20  N 
RECEIVER  OEPTH  -  20  It 
K  •  1.5  It 
k  -  1.0 


Figure  6.  Relative  average  intensity  vs  range  for  a  fixed  depth  using  a  single  source  and 
point  receiver. 
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SECTION  11.  SUMMATION  OF  PRESSURE  WAVES  FROM  MULTIPLE  SOURCES 


Some  of  the  theory  for  the  single  source  problem  is  not  directly  applicable  or  is  not 
practical  for  the  solution  of  a  multiple  source  problem,  particularly  if  the  solution  is  to  be 
obtained  on  a  computer  for  a  large  number  of  sources. 

In  equation  (6)  the  sum  of  two  pressure  waves  was  squared  and  by  the  use  of  trigono¬ 
metric  identities  was  reduced  to  something  which  could  be  time  averaged.  If  the  problem 
required  two  sources,  then  it  would  be  necessary  to  square  the  sum  of  four  pressure  waves, 
to  attempt  to  use  trigonometric  identities  for  reduction,  and,  finally,  to  time  average.  This 
process  is  impossible  when  considering  1000  or  more  sources  in  a  real  sonar  system. 

In  contrast  to  the  trigonometric  solution,  complex  notation  provides  direct  and 
usable  quantities  for  computer  processing.  From  figure  4  it  is  apparent  that  we  may  sum 
the  components  of  as  many  pressure  vectors  as  we  wish.  Simple  squaring  of  the  summed 
components  will  provide  a  squared  pressure  vector  as  shown  in  figure  7. 

The  sum  of  the  pressure  vectors  would  be 

- ► 

(P,  +  P-,  +  P3  +  P4)  =  (A  +  B  +  C  +  D)  +  (E  +  F  +  G  +  H)  i .  (35) 


i 


Figure  7.  Plot  illustrating  the  addition  of  vectors  representing  complex  numbers. 
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To  find  the  resultant  pressure  squared,  simply  square  each  component  and  add. 

(P,  +  Pt  +  P3  +  P4 )“  =  ( A  +  B  +  C  +  D)2  +  (E  +  F+  G  +  H)2.  (36) 

The  angles  of  the  components  A.  B,  C,  etc,  are  defined  by  the  quantity  within  the  brackets 
in  the  pressure  wave  equations  (17-20)  (tuT  -  2rrr/X).  However,  the  time  part  was  factored 
out  in  equation  (20)  so  the  angles  become  (-  2rrr/X).‘  The  time  averaged  relative  intensity 
of  a  number  of  sources  with  direct  and  reflected  paths  can  generally  be  written  as  (reference 
equations  (26)  through  (29)): 


Aj,  Ar.  Bj,  Br .  .  .  will  be  determined  by  the  cardioid  of  the  source  and  the  direction 
to  the  receiver  or  reflection  point,  and  r^j,  r^r,  rgj,  Tgr  ...  are  the  path  lengths  and  will  be 
determined  by  the  geometry  of  the  problem.  Also,  the  minus  signs  in  the  phase  were  left  out 
because  cos(-  0)  =  cos  d  and  sin(-  6)  =  -  sin  0.  The  quantity  is  squared  so  the  sign  has  no 
effect. 


SECTION  III.  SOURCE  CARDIOID 

A  source  may  be  defined  by  a  cardioid  in  polar  coordinates  (A.  vl'  >  as  shown  in 
figure  8.  The  cardioid  is  symmetrical  about  the  boresight  axis.  The  pressure  A  at  a  given 
point  on  the  cardioid  is  a  function  of  the  angle  i//  from  the  boresight  axis  and  the  reference 
pressure  A|  on  the  axis.  Thus, 

(3, 

At  this  time,  the  author  is  unsure  of  the  effects  upon  the  model  of  a  choice  of  one 
source  cardioid  over  another.  Evidence  by  D  Smith  at  ARL/UT  and  others  shows  that  the 
baffle  of  the  transducer  can  be  included  into  the  model  by  a  correct  choice  of  source 
cardioids.  Theory  for  deriving  a  proper  cardioid  is  treated  in  references  3-5.  Reference  5 
uses  an  equation  similar  to 

sin  (kC  sin  i//) 

Ai  -  ~  :  ;  .  13 

1  kfi  sin 


‘Plus  an  additional  “Initial”  phase  for  beam  forming  if  required. 

3.  Summary  Tech  Report  of  Div  6.NDRC,  vol  16,  Scanning  Sonar  Systems,  1946,  chap  9 

4.  Directionality  Patterns  for  Acoustic  Radiation  from  a  Source  on  a  Rigid  Cylinder,  The  Journal  of  the 
Acoustical  Society  of  America,  vol  24.  number  1 .  January  1952,  Laird  and  Cohen 

5.  Technical  Memorandum:  Some  Redundancy  Effects  on  AN/SQS-26  Performance,  Tracor  Document 
#63-233-C.  9  September  1963 
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Figure  8.  Cardioid  in  polar  coordinates. 


where 

k  =  wave  number  (^) 

A 

£  =  element  half  width  . 


Our  model  is  currently  (19  Jan  79)  using 

A  =  s,n  (kism  ^  \1  cos  * )  2  .  (ref  3  and  D  Smith  ARL/UT) 

k£  sin  4/  2 


Equation  (40)  apparently  more  closely  represents  the  form  of  a  beam  generated  from 
a  two-dimensional  element  mounted  on  a  solid  surface  than  does  equation  (38).  In  the  case  of 
a  circular  element  of  10-cm  diameter,  £  would  be  5  cm.  In  the  case  of  a  rectangular  element, 

£  would  be  approximately  half  of  the  diagonal  dimension.  Squaring  the  (1  +  cos  \l/)/2  term 
has  the  effect  of  narrowing  the  cardioid.  Figures  9  through  1 3  show  examples  of  cardioids 
with  £/X  ratios  varying  from  0.25  to  5  based  upon  equation  (40). 


SECTION  IV.  THREE-DIMENSIONAL  DETERMINATION  OF 
PATH  LENGTH  AND  SOURCE  INTENSITY 


GENERAL 


A  Cartesian  coordinate  system  will  be  used  to  define  the  location  of  source  and 
receiver  elements.  Horizontal  and  vertical  angles  are  defined  in  figure  14. 

Sonar  systems  normally  employ  clockwise  notation,  however.  Therefore,  input  and 
output  of  data  will  be  in  the  clockwise  notation  and  converted  to  counterclockwise  notation 
for  the  internal  computational  work  in  the  computer. 

The  point  source  cardioid  of  figure  8  together  with  a  point  receiver  of  omnidirec¬ 
tionality  will  be  placed  in  the  Cartesian  coordinate  system  as  shown  in  figure  15. 
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Figure  14.  Three-dimensional  right-hand  rectangular  coordinate  system  used 
definition  of  horizontal  and  vertical  angles. 


SURFACE 


RECEIVER 


BORESIGHT 


Figure  15.  Configuration  of  source  and  receiver  in  general  case 


The  following  assumptions  apply: 

1 .  The  surface  is  parallel  to  the  x-y  plane  and  can  therefore  be  defined  by  a  single  z 
coordinate. 

2.  The  bottom  (not  shown  in  figure  1 5)  is  parallel  to  the  x-y  plane  as  well,  so  it  can 
be  defined  by  a  single  z  coordinate. 

3.  The  source  cardioid  will  be  in  polar  coordinates  (A,  \p). 

The  following  data  will  be  given  as  input: 

1.  X,  Y,  Z  coordinates  of  the  source 

2.  0,  the  angle  of  the  boresight  (horizontal  plane) 

3.  <p,  the  angle  of  the  boresight  (vertical  plane) 

4.  X,  the  wavelength  of  interest 

5.  C,  the  element  half  width 

6.  Xr,  Y^,  Zj^  coordinates  of  the  receiver 
7. 1  coordinate  of  the  surface 

The  following  intermediate  information  must  be  computed: 

1.  angle  from  the  boresight  to  the  direct  path  between  source  and  receiver 

2.  Ad,  pressure  level  in  the  direction  of  the  receiver  (equation  (40)) 

3.  rd,  distance  from  the  source  to  the  receiver 

4.  1 1>V  Ap  rf,  the  same  computations  as  1  through  3  above  for  the  reflected  path 
(not  shown  in  figure  15) 

Repeat  1  through  4  for  each  source  given. 

The  above  intermediate  data  will  then  be  used  in  equation  (37)  to  compute  the 
resultant  intensity.  By  changing  Xr,  Zr,  or  Z w  the  entire  process  may  be  repeated  to  create 
beam  patterns. 

DIRECT  PATH  COMPUTATION  (*d,  Ad,  rd) 

The  direct  path  distance  rd  may  be  easily  calculated  from: 


rd  =  KXr  -  XS)2  +  (Yr  -  Ys)2  +  (ZR  -  ZS)2] 
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If  the  boresight  is  defined  as  vector  A  and  the  path  from  source  to  receiver  is  defined 
as  vector  B,  then  from  vector  theory 


A  •  B  =  |A|  |B|  cos  i//j. 


(42) 


Therefore,  may  be  obtained  by  writing  the  vectors  for  A  and  B.  By  definition 
A  =  ajx  +  asy  +  a3Z. 

If  we  let 


then,  from  figure  1 5  , 
as  =  I  •  tan  0 
and 


tan  <p 
cos  0 


or. 


A  =  x  +  tan  0 y  + 


tan  0 
cos  0 


7.. 


Also,  by  definition. 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


which  may  be  reduced  to 


IAI  =  1 

cos  d  COS  <t> 

(50) 

To  continue: 

B  =  bjx  +  bsy  +  b^z. 

(51) 

Therefore, 

B  =  (xr  -  Xs)x  +  <yr  -  Ys)y  +  (ZR  -  Zs)z 

and 

(52) 

IBI  =yJ~(XR  -  XS)-  +  (Yr  -  YS)-  +  (ZR  -  ZS)2  . 

(53) 

Since 

A  ■  B  =  a  |  b  |  +  aibs  +  3363  , 

(54) 
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then,  from  equation  (42), 


(XR-XS)  cos  0  cos  0  +  (Yj^  —  Y g)  sin  6  cos  0  +  (Zj^  —  Z§)  sin  0 


(XR-XS)"  +  (YR-YS)-  +  (ZR-ZS) 


Finally,  from  equation  (40), 


sin  (kC  sin  \p^)  f  1  +  cos  i 


REFLECTED  PATH  COMPUTATION  (*r,  Ar,  rf) 

The  reflected  path  solution  may  be  obtained  by  referring  to  figure  16 


SURFACE 


RECEIVER 
R  (Xd.Y„,Z„I 


/  SOURCE 

A3 

'  (XS-YS'Zs' 

PARALLEL 
TO  X-Y  PLANE 


BORESIGHT 


Figure  16.  Reflected  path  from  source  to  receiver. 
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The  reflected  path  is  SQR.  The  surface  is  at  Z^.  0r  are  all  equal  for  obvious  geometric 
reasons.  The  reflected  path  distance  is  then  equal  to  SR  ,  or 


rr  =  (SB2  +  R'  B2)‘/2 


SB2  =  (Xr-Xs)2  +  (Yr-Ys)2 


R'  B"=  [ZR  —  Zs)  +  —  ZR)1  , 


rr  -  [(XR  ~  XS)_  +  (YR  “  YS)_  +  (2Zu  "  ZR  “  ZS)2]  • 


0r  is  defined  as  the  angle  between  the  boresight  and  the  reflected  path  in  the  direc¬ 
tion  of  SQ  of  figure  16.  As  in  the  previous  section,  two  vectors  must  be  defined.  Vector  A 
is  the  boresight  vector  and  is  identical  to  that  of  equations  (47)  and  (50)  referenced  to 
figure  15.  Repeated: 

tan  d> 

A  =  x  +  tan  Oy  + - -z  (62) 

cos  0 


cos  0  cos  0 

Vector  B  may  be  similarly  expressed  in  terms  of  0r  and  0T  as  shown  in  figure  16. 

tan  0r 

B  =  x  +  tan  0.  y  + - - — z  . 

r  cos  6r 


From  figure  16, 


„  (YR_YS) 


cos  6  r  = 


XR-XS 


and  tan  0r  - 


[(Xr-xs)2  +  (Yr~ys)“] 

(2Zco~ZS~ZR) 

|(xr  -  XS)2  +  (yr  ~  yS)2J  ' 


->V/2  ’ 


Substituting  into  equation  (64),  we  have 

_  (YR~  YS}  (2Zoj~ZS~ZR) 

=  x  +  (xR-xs)y+  (XR-XS)  z 
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